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Abstract 

Diffractive deep inelastic events with a large rapidity gap are analyzed 
by using a Regge model for the pomeron flux and a gluonic content for the 
pomeron. Contrary to the expectations, the simplest assumption for the 
pomeron trajectory gives the best agreement with the data on the ratio of 
diffractive to the total number of events. In this case the main properties of 
the model are described by an analytic expression. 
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1. Introduction 

Recently events with a large rapidity gap have been revealed at HERA [1] that 
are not expected on the basis of standard deep-inelastic Montecarlo programs. In 
this paper we show that the experimental results are compatible with a theoretical 
model based on the diffractive nature of these events. 

By definition, the differential cross section of a diffractive deep-inelastic process 



where y is the pseudorapidity, x = /2q-p is the fraction of the proton momentum, 
^ = 1 — xp = p'zIPz and t is the squared proton momentum transfer, t = [p—p'Y. 
The 2;-axis points along the direction of the incoming proton. Note that > x and 
that, for a diffractive event | t |< (few hundred MeV^) and ^ 1. We remind also 
that at HERA, {k -\- pY = \/s = 29QGeV and the pseudorapidity of the smallest 
detector angle is rj = 4.3,0 = 1.5°. The cut in rj^ax < 1-5, distinguishes events 
with a large rapidity gap and is equivalent to ^max ^ 0.06. Due to acceptance cuts, 
the maximum value of ^ is = 2.0 X 10~^ for ZEUS [1]. 

Eactorization 



most of the calculations of diffractive deep inelastic scattering. We believe that this 
assumption does not need modifications. The distribution of the invariant mass 
of the hadronic system produced is peaked at small mass values and the triple- 
pomeron vertex, that appears only in the triple- Regge limit, does not contribute in 
this case. 



(Fig. 1) is 







of the structure function, where -Pp/p is the so called pomeron fiux (see below) and 
Gq/p IS the pomeron structure function, is an important assumption [2-5] adopted in 
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Figure 1: Diagram for the diffractive deep inelastic scattering. 

2. The pomeron flux 

The pomeron flux 

'''^ ~ (TtotiPp) didt 
has been introduced [6] in the description of the diffractive dissociation. 

We shall investigate the role of the t-dependence by using also non linear pomeron 
trajectories, fltted earlier [7] to elastic pp and pp scattering. More important, we 
shall test also the role of the energy dependence in the pomeron flux. We remind 
that in earlier papers a constant value for the pomeron-proton total cross section 
was used (e.g., 2.3 mb in [8] or 1 mb in [6]). We shall include the energy dependence 
in the pomeron flux explicitly. 

We start from the following simple form for the pomeron flux -Pp/p, borrowed 



2 



from our experience in describing high energy elastic hadron scattering 

Fp/,(e,t) = [exp(i?a(t))r"(*)]2e, (3) 

where a[t) is the pomeron trajectory. In the case of a unit intercept hnear trajectory 
a[t) = 1 + a't and for B = 10 [9], the expression (3) is numerically equivalent in 
the t-range under consideration to that for the pomeron flux used by Donnachie 
and Landshoff [10], even if the analytic form of the residue functions is different. 
By using a non linear pomeron trajectory, we shall study the effect of the non- 
exponential behaviour in t as well . 

The variation with ^ of the pomeron flux may be studied by letting a(0) to be 
slightly beyond I ("supercritical" pomeron) or by using a dipole pomeron instead 
of Eq.(3), resulting in a logarithmic rise of the pomeron flux in ^. 

Practically, the pomeron flux cannot be completely isolated from non leading 
contributions that modify Eq.(3) as 

FiU) = Fp/,iU) + aRFn/^iU) , (4) 

where 

i^i^/.(e,o = M''"""^*^^r""^*¥e, (5) 

aR[t) is an effective Regge trajectory and is the relative contribution of the 
"Reggeons". In this paper we will not consider the contribution from "secondary 
Reggeons" . As will be shown elsewhere their presence does not change sensibly our 
conclusions. 

3. Evolution of the pomeron structure function 

To calculate the effect of the Q^-evolution of the pomeron structure function 
we use a very simple input distribution corresponding to a pomeron made of a 

3 



small number of ocean gluons only. In this case, and for non-large values of 
(Q^ = 5 lOGeV^)^ we assume a gluon distribution of the form [8] 

zG{z,Ql) = 2{l-z) , (6) 

satisfying the momentum sum rule 

C dz zG{z,Ql) = l . (7) 

JO 

We notice that Eq.(6) does not contrast with the experimental results [11] since, 
till now, the data cannot discriminate between the two choices 1 — z an zil — z). 
Radiative corrections to the pomeron gluonic ladder suggest that the form (6) could 
be preferable but, at this stage, it can be considered as an ansatz. 

The solution of the evolution equation in a wide range of z necessitates cumber- 
some calculations and assumptions that make it difficult to get further transparent 
results for the diffractive DIS. In particular, it becomes problematic to take into 
account gluon recombination effects for the initial distribution (6) [3]. Expression 
(6) is fairly simple to enable an estimate of the evolution that will reproduce the 
main features of the gluon distribution in the kinematical region of interest, where 
z = J assumes not too small values [z > 0.01), being x ~ 0.001 and !^max ~ 0.06. 

Thus, we start from the distribution 

zG{z,Q') = a{Q'){l-zf(^') , (8) 

in the kinematical region under consideration, with a{Q^) = 2 and h{Q^) = 1. 

The expression for b[Q^) can be easily obtained near 2; = 1, when quarks are 
neglected, by going to the moments. For z ^ 1 the threshold behaviour of zG{z^ Q^) 
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is implemented by the large n dependence of the anomalous dimension at the one- 
loop approximation [12]. We get 

12 

HQ') = 1 + giv, ^ . (9) 

3 

with 



s = ln{- 



Ql = ^GeV^, A = 0.2G'ey. 

Increasing from Q\ will diminish the value of the momentum integral (7). If we 
write = 2exp(as) the bound 

< 1 + \_^ s (10) 

3 

must be satisfied. 

According to Ref. [13] the momentum integral is reduced by 10% at = 
50 GeV^ for the distribution we are considering. This condition determines a and 
the approximation we will use in this paper. With three fiavors, Nj = 3, we obtain 

HQ') ~ 1 + f , <Q') ~ 2e^^' (11) 

and the resulting evolution for z > 0.01 is shown in Fig. 2 Let us now calculate the 
probability that a pomeron contains a parton with a momentum fraction z from 
the evolution equation 

dq{z, Q"^) a,s /-i dy 



ain(gvA') 2^-^- y 

with 



^ / -^yGiy,Q')P,,iz/y) , (12) 



PUz/y) = 1/2 -z/y + z^y' 
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Figure 2: z -dependence of the gluon distributions: the initial distribution at = 
5 GeV^ (dotted line) is evolved to = 50 GeV^ (dashed line); a sketch of th 
Ingelman-Prytz solution [2] is also shown (dotted-dashed line) together with th 
standard DGLAP evolution (full line) evaluated in Ref. [2]. 
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and 

yG{y,Q') = a{Q'){l-yf(^') , 
where a{Q^) and b[Q^) are known. 

Within the approximation suggested by Ingelman and Prytz [2], we obtain from 
Eq.(12) for Nj = 3 (see Appendix) the crude estimate 



|(i-.)^+W) 



■ 25-1 (2-5X6-1)^, 
. 6 6 



^b (. , .^ (&-l)(2-&) 2 



l + {h-l)z-'- '-^ '-z']2Fi{l,h+l-h + 2-l- z 



(13) 



2(6+l)V ' ' 3 
Other more rigorous approaches exist, see for example Ref. [3], but at the price of 

introducing one new parameter. 
4. A simple model 

Before starting with the numerical calculation, it may be useful to understand 
better the problem in a simplified model. Eq.(13) shows that a simple input for the 
pomeron fiux will allow for an analytical solution. In Eq.(3)) we introduce a linear 
trajectory with a' = 0.25 GeV~^ and write 

Fp/,(e,t) = c(e^(i+«'*)r"^*))^e = 

normalized according to Ref. [6] at t = to 

Fp/p(^, 0) = 3.536/^ ^ Ce^^ = c = 3.536^6^-^ , (15) 

whence 

FpMe,t) = p«'("«)*. 



Integration over t gives 



Now, since the numerical value of the exponent 

2a' {B - In 

is typically 0.4 for to ~ 0.05GeV^,a' = 0.25^6^-2 and ^ = 0.0009 (note that the 
integration over ^ extends from x to ^ 0.02 according to the acceptance cut in 
ZEUS [1]), we can approximate Eq.(16) as 

r Fp/,(e,tMt^^. (17) 

J-to ^ 

To evaluate the relative contribution of diffractive deep inelastic events at HERA, 
we first have to calculate the quantity 

•J —Iq J X J X t, 

8 . ,^2> J „ .Si-l {'i-b){b-l) 



xi(,o \32; 6 6 



^ , .^ (&-l)(2-&) 2 



2(^^l)(l + (&-l>-' i '-z\F,{l,h+l-h + 2-l-z)j . (18) 

Performing the integration in the right hand side of Eq.(18) (see Appendix for this 
purpose) leads to the result 



8 / \b{Q^)+2 

L{x,Q')c:,-ctoa{Q')[l-x/Co) 
{h{Q^) - l/2)(6(g2)2 _ 1/2) (6(g2) _ ^^^^^q2^ _ 3) 



ln(^o/a:) 

(1 - ^/eo) 



3(6(g2) + 2) 6 
KQ") ,eo, , {h{Q')-l){h{Q^)-2) 



2{l-xl(,)\' Co ^x'^ 3(6(g2) + l) ; 



(19) 
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where c = 3.536^6^-2,^0 = 0.02, a(g2) = 2e'/^,b{Q'^) = 1 + 4s/3. By using the 
parametrization (10) of Ref. [14] for ^2(2;, Q^) we are now able to calculate the ratio 

r - (20) 

The results of our calculations for various values of x and for in the interval 
between 10 and 100 GeV^ and to = 0.05 GeV^ are shown with full lines in Fig. 3, 
where also the ZEUS data [1] are shown. Let us notice that we can arrive to the 
same results normalizing Eq.(14) according to Berger et al. [8] and choosing to ~ 0.1. 
The independence of our finding from the parametrization chosen is exhibited in 
Fig. 4 where, by using Fq.(19), we fit preliminary data for the F2^^^^'{x^ Q^) [15]. 

5. Alternative models 

A numerical analysis has been performed with the non linear trajectory 

a(t) = ao + at - 7(yt7^ - Vtl) (21) 

where 

a' = 0.25 GeV-^ , 

7 = 0.03 GeV-^ , 

t, = Ami - 0-078 GeV"^ . 

When t << tg the trajectory (21) reduces to a linear one but, beyond the threshold 
ts, it describes the transition to the Orear regime of the amplitude [16]. The virtue 
of this parametrization is that it accounts also for the small-t structure observed 
in the diffraction cone in high energy hadron scattering (the so-called "break"). In 
ref. [16] more details and earlier references are given. 
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Figure 3: The ratio r versus for x = 0.0018 (a), 0.0009 (b), 0.0045 (c) and for 
to = 0.05 and = 0.02, compared with the ZEUS data [1]: linear model (full lines) 
and supercritical pomeron with non linear trajectory (dashed lines). 
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Figure 4: F^'^ ^'■{x.Q^) versus log^o{x) for = 8.5 GeV^ (a), = 15 GeV^ (b), 
= 30 GeV^ (c) and = 60 GeV^ (d), as given in Eq.( 19) of the text, conapared 
with the data from Ref. [15] 
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In a first instance tlie value of ao lias been set equal to one. The calculation 
provides also a test for the simple model of Sec. 4. The difference in the results as 
compared to the case of linear trajectories was expected to be small, since the |t| 
integration range was restricted to a limited region near the forward direction and 
also because the non linear effects in the trajectory in this region are small anyway 
(numerically, 7 <^ a'). Surprisingly, the calculated effect in r is considerably large 
and, in general, it tends to deteriorate the agreement with the available data. 

The same is true for the pomeron intercept. Since the integration range in the 
energy variable (pomeron flux) is limited, one would expect that the choice of the 
pomeron singularity (i.e. the value of the pomeron intercept) has little effect on 
the final results. In fact, this is not the case: the introduction of a new parameter, 
5 = a(0) — 1 (in our calculations we put 5 = 0.1) again deteriorates the results. 
Finally the results are not better in the case of a "mixed" option, i.e. when both 
the trajectory is non linear and the pomeron intercept is "supercritical". As an 
example, we show in Fig. 3 (dashed lines) the result of the computer calculation 
for the ratio r with a supercritical pomeron, 5 = 0.1 and the non-linear trajectory 
(21). 

It should be stressed that by working with non linear trajectories and a su- 
percritical pomeron, we have not introduced free parameters - their values were 
taken from earlier fits of elastic hadron scattering [7] and were not varied here. 
Since, however, those fits and the relevant parameters themselves may vary, this 
point needs further clarification. Here we can only assert that good agreement has 
been achieved in the "minimal" model of a linear trajectory and simple pomeron 
pole with a unit intercept (constant asymptotically total cross sections) and that 
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the final results are sensitive even to small variations of the parametrization of the 
"pomeron flux" . 

6. Conclusions 

The aim of the present paper was to study the pomeron internal structure by us- 
ing diffractive deep inelastic scattering date from HERA [1]. We have constructed a 
model for the diffractive deep inelastic scattering based on the formalism developed 
earlier [2] as well on our present knowledge about the elastic scattering amplitude 
(pomeron exchange) and the assumed pomeron structure (distribution functions). 
By comparing the calculated dependence (at various fixed values of x) of the 
relative contribution of the diffractive DIS events, Eqs.(17) and (18) to the HERA 
data, the role of the various components of the input, namely the pomeron structure 
function, the form of the pomeron trajectory and its intercept, as well as the effect 
of various approximations used in the calculations were tested. 

We have tried different inputs as regards the form and the intercept of the 
pomeron trajectory. A non linear trajectory has been also considered that, for 
hadronic reactions, describes the small-t "break" in the differential cross-section. 
Besides the unit intercept, the case of a supercritical one, with a(0) = 1.1 has been 
analyzed too. 

While the results of the calculations for different inputs, including also mixed 
options, sensibly differ, the main conclusions are quite general. We found noticeable 
dependence on the pomeron intercept, which means that, in spite of the limited 
interval for the energy (pomeron flux) variation, its rate has important consequences 
for the ratio r. 
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The introduction of non-linear trajectories does not improve the resulting be- 
haviour of the ratio r. A better choice of the parameters may improve the agreement 
with the data, but anyway this is only a small effect as compared to the available 
freedom in choosing the pomeron structure function. 

While the variation of the pomeron flux is meant mainly "to fix normalization" , 
the objective of the present paper - like in similar investigations by other authors [2, 
3, 6, 8, 10] - was the right choice of the pomeron structure function, i.e. the 
deduction of the pomeron internal structure from the data. We have assumed that 
the pomeron is made of a small number of gluons distributed according to Eq.(6), 
with the parameters a[Q^) and b[Q^)^ calculated from the evolution equation, given 
by 

6(g2)~i + 4^, a(g2)~2et 

and normalized according to a((5o) — 2, &(Qo) — 1- Note that the above distribution 
is close to the one calculated by G.Ingelman and K.Prytz [2] taking into account 
gluon recombination. Any modification of the pomeron distribution by using e.g. 

zG{z,Ql) = Q{l-zY 

(corresponding to the presence of many gluons inside the pomeron) instead of Eq.(6) 
and its QCD evolution, sensibly affects the resulting behaviour of the ratio r. The 
large errors in the experimental diffractive structure functions do not allow, however, 
to exclude this possibility. The agreement with the data of Ref. [1] of the results 
of our calculations is quite satisfactory since, for a given choice of the pomeron 
structure function, practically no parameters were introduced (those used in our 
calculations were fixed elsewhere - for example the pomeron trajectory was adjusted 
earlier to fit hadronic cross sections). This agreement can be made even better by 
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letting some of the free parameters to vary and to be fitted to the diffractive deep 
inelastic data, which however are still quite preliminary and subject to further 
improvement. 

To summarize, we consider our model as a laboratory for studying the pomeron 
structure function by comparing the calculated results with the experimental data. 
In this paper we made the first step by proving the stability of our model against the 
variations of its component and tentatively using a rather simple gluon distribution 
function for the pomeron internal structure. 

Acknowledgement: We thank M. Arneodo, A. Garfagnini, L. lannotti and A. 
Solano for useful discussions on the subject of this work. One of us (L. L. J.) 
is grateful to the Dipartimento di Fisica dell'Universita della Calabria and to the 
Istituto Nazionale di Fisica Nucleare - Gruppo coUegato di Cosenza for their warm 
hospitality and financial support. 

6. Appendix 

Let us start from Fq.(ll). By using the Ingelman-Prytz approximation [2], for 
Nf = 3 we get 

Ztt Jz y 

2/9a{Q') f dy/y'{l - yf^'Hl/^ - z/y + z'/y') . 

J z 

The integration can be done explicitly yielding 

/ = dyly^l - y)^(^')(l/2 - z/y + z'/y') = 1/21, - zh + z'h , 

J z 
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where 

By applying to 2Fi{n^ b -\- l;b -\- 2;1 — z) the recurrence formula 

F{n) = -^—{{b + 1)^-"+^ + (n - 6 - 2)F(n - 1)) , 
n — 1 



we get 



^ ^ ('(&+!) ^ (&+l)(2&-l) _ {2-b){P-l)z 
3z Q 6 



-1/2(6 + 6(6 - l)z - 6(6 - 1)(2 - 6)^73)2^^1(1, 6 + 1; 6 + 2; 1 - z)^ 



(1-^)^+1 



6+1 

The integration in z may be easily done by using the form 



,Fi(l, 6 + 1; 6 + 2; 1 - z) = (6 + 1) /' dt- 

Jo 1 



t{l-z) 

We can approximate the result with the limit 



lim$(z,l,t;)/(-ln(l - z)) = 1 



where 



To get 



^z,l,v) = J2{v + n)-'z'' . 



n=0 



f 

it is sufficient to multiply zq[z^ Q^) by 4/3, {N j = 3). 
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